Quasimodes of a chaotic elastic cavity with increasing local losses 
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We report non-invasive measurements of the complex field of elastic quasimodes of a silicon wafer with 
chaotic shape. The amplitude and phase spatial distribution of the flexural modes are directly obtained by Fourier 
transform of time measurements. We investigate the crossover from real mode to complex-valued quasimode, 
when absorption is progressively increased on one edge of the wafer. The complexness parameter, which char- 
acterizes the degree to which a resonance state is complex- valued, is measured for non-overlapping resonances 
and is found to be proportional to the non-homogeneous contribution to the line broadening of the resonance. A 
simple two-level model based on the effective Hamiltonian formalism supports our experimental results. 
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A closed quantum system is fully defined by a set of 
eigenenergies and orthogonal discrete states. When the sys- 
tem is coupled to the environment, e.g. through leakage at the 
boundaries, mode lifetime becomes finite. Consequently, the 
spectral widths of the resonances broaden and are no longer 
isolated. While the statistical properties of the spectral widths 
of chaotic wave systems was systematically analyzed in the 
regime of isolated resonances 01 Hi, in the overlapping regime 
UIU] and in the strong overlapping regime 0], the effects of 
the leakage on the wave function statistics remain an issue 
|@, 0]. For a wave system whose closed limit displays time re- 
versal symmetry, the eigenfunctions become complex- valued 
quasimodes due to the presence of currents, the standing- wave 
component being progressively replaced by a component trav- 



. Such a complexness 
10] and long-range cor- 
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eling toward the system boundaries 
may reveal itself in current density 
relations of wave function intensity 
non-real nature of the field appears also in various domains 
of wave physics. It is an essential ingredient in the theory 
of lasing modes, which induces an enhancement of the line 
width as pointed out by Petermann lfl3ll and studied in de- 
tails for chaotic lasing cavities fl4ll . The complex nature of 
the field was also discussed recently in the context of disor- 
dered open media fl5ll and in diffusive random lasers fl6h . To 
quantify the complexness of the field, it is convenient to intro- 
duce the ratio of the variances of the imaginary and real parts 
of the field as a single parameter: the complexness parame- 
ter q 2 f\j\\ . Measuring q 2 requires the complete knowledge 
of the spatial distribution of the field, including its amplitude 
and phase. Indirect estimation of the complexness parameter 
assuming ergodicity over several resonances was reported in 
But, to the best of our knowledge, direct measurements 
for a given quasimode has not been reported. 

In this work, we measure non-invasively Lamb waves prop- 
agating on a doubly-truncated silicon wafer and measure the 
effect of increasing losses on spectral and spatial characteris- 
tics of the modes. Leakage channels are progressively opened 
by sticking absorbent strips with increasing dimensions along 
one edge of the chaotic sample. For each configuration, the 



acoustic field spatial distribution, the phase probability distri- 
bution, the complexness parameter q 2 and the spectral width 
T are measured for individual non-overlapping resonances. 
Experiments show a simple proportionality between q and T 
when the losses are varied. Experimental results are found to 
be consistent with a simple analytical 2-level model based on 
the scattering approach of open systems. 

The complexness parameter of the nth quasimode \I/ n is 
defined by: 



(Im(*n) 2 ) 
(Re(^n) 2 ) 



(1) 



where the triangular brackets denote the spatial average. It 
is worth noting that q 2 is equal to for a closed cavity with 
no currents and tends to 1 for a pure traveling-wave in open 
space. 

The derivation of the probability distributions of the com- 
plexness parameter in the perturbative regime was derived in 
fl9l1 using the effective Hamiltonian formalism and applying a 
random matrix approach to open systems (see J20l UM for re- 
views). While a TV-level model with N — > oo is needed for the 
probability distribution, a 2-level model |22l [23|] is here suffi- 
cient to consider the relationship between the spectral widths 
and the complexness parameter. We start from the effective 
Hamiltonian H e ff = H — iVV T /2, where H is the Hamil- 
tonian of the closed system modeled by a 2 x 2 matrix and 
iVV T /2 is an imaginary potential describing the coupling to 
the environment in terms of M open channels. The 2 x M 
matrix V contains the coupling amplitudes V£ which couple 
the nth level to the cth open channel. As a result of the non- 
hermiticity of the effective Hamiltonian, its eigenvalues and 
eigenvectors are complex. The eigenenergies of H e f f reads 
£1,2 = ci,2 — f Ti,2, where ei 5 2 and Ti^ are, respectively, the 
two eigenenergies and the two spectral widths of the 2-level 
model. In the eigenbasis of H the effective Hamiltonian is 
written as 



H eff = 









E 2 



1 An Ti 2 
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(2) 



where Ei^ are the eigenenergies of H (E 2 > E\ is assumed) 



and T np = Y^=i ^n^p- As we f° cus on the isolated res- 
onance regime, the imaginary potential may be viewed as a 
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perturbation of the Hermitian part and the effective Hamilto- 
nian © can be easily diagonalized through a first order pertur- 
bation theory ll33ll . One obtains straightforwardly the eigenen- 
ergies = £1,2, the spectral widths = X^i^i^) 2 ' 
while the perturbed eigenvectors and ^2), read, in the 
basis {|1>,|2» of H: |^) = |1) - if\2) and |^ 2 > = 
|2)+i/|l),with/ = r 2 i/(2(E 2 - J Bi)) =T 12 /(2(E 2 -E 1 )) 
(in this model Ti2 = r 2 i). 

In the 2-level model, q 2 is obtained using (Q} by replacing 
the spatial average by an average over the components of the 
eigenvectors and 1^). It reads q 2 = r 21 /(4(£^ 2 — Ei) 2 ), 
for both resonances . For a uniform increase of the inhomo- 
geneous losses [24], each coupling amplitude increases in the 
same way: V — > y/vV, where y/v characterizes the enhance- 
ment of the losses. As a result, the spectral width and the com- 
plexness parameter depend on v as: r — > vT and q 2 —> v 2 q 2 
implying a linear relation between T and q for a given mode 

9 " E c (Kf) 2 2(E 2 -E 1 )- {) 

The relationship can be written under the general form q n = 
P n T n , where the slope f3 n depends on the resonance because 
of the fluctuations of both coupling amplitudes and spectrum 
i25ll . In the limiting case of a large number of weakly cou- 
pled channels (M -> 00, a 2 -> and (r) = Ma 2 fixed) the 
term depending on the coupling amplitudes is proportional to 
l/y/M as reported in W22I1. such that the fluctuations of the 
slope are only due to the spectrum. In the following we will 
only focus on the proportionality between q n and Y n consid- 
ering a 2D-chaotic acoustic cavity. 

Two kinds of elastic waves can propagate in thin plates: 
horizontally polarized shear waves (SH) and Lamb waves, 
which are a combination of vertically polarized shear waves 
(SV) and longitudinal waves (L). SV- and L- waves cou- 
ple with each other on plate/air interfaces and provide sym- 
metric and antisymmetric displacements i26ll . The only 
guided modes which exist at all frequencies are the zero- 
order symmetric #0 and antisymmetric Ao modes, also called 
extensional and flexural modes in the low frequency limit 
(wavelength^>plate thickness). In the low frequency range, 
Ao is highly dispersive. Here we consider Lamb waves prop- 
agating on 2" silicon wafers of thickness h = 380/im. The 
initial silicon wafer is cut at R/2 to form a D-shape plate. 
The associated classical billiard is known to display chaotic 
dynamics of ray trajectories i27ll . This chaotic shape ensures 
mode ergodicity as opposed to integrable cavities, such as cir- 
cular ones, where modes have regular patterns. The modal 
statistics in such a chaotic cavity is essentially Gaussian with 
spectral repulsion between nearest eigenfrequencies B28I1 . An- 
other cut perpendicular to the first one is made at \ /r 3R/2 to 
break the remaining symmetry. The resulting shape is shown 
in inset in Fig. [4] 

The acoustic waves are optically excited using a frequency- 
doubled Nd:YAG Q- switched laser (Quantel Ultra) operating 
at wavelength 532 nm and 8 ns pulse width. The laser pulse 
heats the surface of the sample and creates a few /is-long 
acoustic pulse by thermoelastic effect. The beam is focused 
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FIG. 1 : Time measurement at one point of the sample. Inset: (top) 
beginning of the signal, (bottom) corresponding spectrum. 



down to a 200 /im-diameter spot, which corresponds to a flu- 
ence of 0.48 J/cm 2 . This is below ablation threshold 112911 and 
is a good compromise between bandwidth and signal-to-noise 
ratio. 

On the other side of the wafer, a Thales SH-140 interfer- 
ometric heterodyne optical probe measures a time response 
proportional to the plate normal displacement at one point 
over a wide bandwidth [20 kHz-45 MHz] . The resolution of 
10 -4 A/\/Hz allows a sensitivity to displacements of the or- 
der of a few angstroms L30]. As our probe is only sensitive 
to the normal displacement of the plate, we detect preferen- 
tially Aq. The time sequence of 500 thousand data points is 
recorded by a 400 MHz-bandwidth Lecroy WaveSurfer digi- 
tal oscilloscope triggered by the laser pulse at a repetition rate 
of 20 Hz, and averaged over 100 traces. The laser excitation 
and the optical detection provide a totally non-invasive setup 
to investigate the dynamics of free acoustic wave propagation. 

The sample is supported horizontally by three pins. We 
checked that their impact on acoustic propagation is negli- 
gible. The main source of losses comes from coupling with 
air. To scan the entire wafer, a XY-translation stage is used 
with a 500 /im-spatial resolution, while the optical probe re- 
mains still. Mirrors are attached to the stages to ensure that 
the excitation hits the sample at the same position after each 
displacement. Figure [T] shows a typical time-signal recorded 
at one point of the sample. The dispersive predominant anti- 
symmetric Ao mode reaches the optical detector after a few [is 
and is subsequently reflected at the edges of the wafer. Nei- 
ther the symmetric So Lamb mode foregoing Ao, nor the SH- 
mode are detected by the optical probe. The entire multiply- 
scattered exponentially decaying signal lasts several millisec- 
onds. Time records are stored for each value of X and Y. 
Hence we reconstitute the space-time map of the normal dis- 
placement field at the surface of the wafer. A speckle-like 
spatial distribution is rapidly reached as a result of the chaotic 
geometry of the doubly-truncated wafer 13 ill . 

Each time-record associated with each point on the wafer is 
Fourier-transformed to obtain the spatial distributions of the 
real and imaginary parts of the acoustic field as a function of 
frequency. An example of a power spectrum measured at one 
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FIG. 2: (Color online.) Real and imaginary component of a quasi- 
mode at its central frequency for: (1) no absorber, (2) 1 mm- thick, 
45 mm-long and 1 mm-lage strip, (3) 2 mm-thick, 45 mm-long and 
1 mm-broad strip, (4) same as (3) with a 2 mm-broad strip, (5) same 
as (3) with a 3 mm-broad strip. 



point is shown in lower inset of Fig. [T] Each peak corresponds 
to a vibration mode of the plate. The intensity of each peak 
depends on the overlap of the mode with the excitation and 
detection. For a well isolated resonance, one obtains the maps 
of the real and imaginary parts of the corresponding mode. An 
example is shown in Fig. [2] We measure the complexness pa- 
rameter and the spectral width for individual non-overlapping 
modes in a set of experiments where absorption is increased 
locally on one edge of the sample. Here we are interested in 
relating q to T, which comes down to comparing the change in 
the spatial nature of the mode with its spectral characteristics. 
This is accomplished by sticking absorbent strips with differ- 
ent dimensions on the wafer's longest edge, as shown in inset 
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FIG. 3: (Color online.) Phase probability distribution of the quasi- 
mode for each of the 5 configurations described in Fig. [2 compared 
with Eq. $4^ with values of q 2 obtained from the ratio of the vari- 
ances of the imaginary and real parts of the quasimode at its central 
frequency. (Bottom right) Corresponding spectra. 



of Fig. SI The absorbent is a PolyDiMethylSiloxane (PDMS) 
elastomer (Dow Corning, Sylgard 184). A 10:1 PDMSxross- 
linker mixture was stirred and degassed then cured at 80 °C 
for 2 hours. While the width and the thickness of the polymer 
strips are varied, the length is kept constant to fix the num- 
ber of opened channels. The degree of acoustic absorption 
of each strip is not fully controlled and depends for instance 
on the way the strip is stuck on the wafer. However, the rel- 
evant parameters are measured for a given value of the ab- 
sorption whatever it can be. For each configuration, the spec- 
trum around 47 kHz where an isolated resonance has been 
identified is plotted in Fig. [3j As the absorption increases, 
the central frequency of the resonance is shifted toward lower 
frequencies while its spectral width broadens. The real and 
imaginary parts of the quasimode are represented in Fig. [2] for 
each sample configuration. The increasing absorption results 
in a global damping of the real part, whereas the imaginary 
contribution increases on average. This is accompanied by a 
progressive deformation of the field spatial distribution. 

To obtain the complexness parameter, modes with Gaus- 
sian statistics are selected, excluding for instance scarred 
modes B32I1 . and a phase rotation is performed on the mea- 
sured field ^ n : \I/ n = e ia "^) n in order for the real and imag- 
inary parts to be independent variables The phase a is 
unique and is fixed by the constraint (Im(\I/ n )Re(\I/ n )) = 
resulting in a probability distribution of the phase cp — 
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arctan(Im(^)/Re(^)) given by O: 



1 



27T COS 2 if ■ 



sin 2 (/? 



(4) 



This expression is peaked around and 7r for purely standing 
waves in a closed cavity. As channels are increasingly opened, 
the phase probability distribution broadens as seen in Fig. [3j 
corresponding to a growing traveling-wave component of the 
mode. The linewidths are obtained by fitting the spectra with 
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FIG. 4: Square root of complexness parameter versus spectral width 
measured for the quasimode at its central frequency for the 5 config- 
urations described in Fig. [2] and linear fit to the data points. Inset: 
doubly-truncated 2"-diameter silicon wafer. 



a complex Lorentzian L(u) — C / (uo — uoq + iT), where the 
central frequency c^o, the spectral width F and the complex 
constant C are fitting parameters. The complexness parame- 
ter is obtained from the ratio of the variances of the imaginary 
and real parts of the quasimode at its central frequency. Figure 
01 shows a linear dependence of q vs. T. The deviation from 
the linear fit q = 5.54 x 10" 4 - 0.04T is below 3%. The com- 
plexness parameter extrapolates to for r=70.5Hz. This value 
corresponds to the homogeneous contribution to the losses re- 
sulting from uniform surface-coupling with air. This losses 
contribute only to a uniform decay rate of the mode. Only 
losses localized at the sample edges contribute to the complex 
nature of the mode. Different slopes are found for different 
modes. A systematic study over a larger number of quasi- 
modes should allow a statistical exploration of the spatial and 
spectral characteristics of the resonance states and should pro- 
vide for an experimental check of the statistical features of the 
complexness parameter lfl9ll . 
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